Preprint typeset in JHEP style - HYPER VERSION 



hep-ph/1001.4878 



Energy-independent new physics in the flavour ratios 
of high-energy astrophysical neutrinos 



M. Bustamante"' 6 , A.M. Gago", C. Pena-Garay c 

a Section Fisica, Departamento de Ciencias, Pontificia Universidad Catolica del Peru, 
Apartado 1761, Lima, Peru 

b Theoretical Physics Department, Fermi National Accelerator Laboratory, Batavia, LL 
60510, USA 

c Lnstituto de Fisica Corpuscular (LFLC), Centro Mixto CSLC-UVEG Edificio 
Lnvestigacion Paterna, Apartado 22085, 46071 Valencia, Spain 



E-mail: nbustamante@pucp . edu .pe , agago@pucp . edu.pe , carlos . penyaOif ic .uv.es 



Abstract: We have studied the consequences of breaking the CPT symmetry in the 
neutrino sector, using the expected high-energy neutrino flux from distant cosmological 
sources such as active galaxies. For this purpose we have assumed three different hypotheses 
for the neutrino production model, characterised by the flavour fluxes at production <fP e : 
4> : <\P r = 1 : 2 : 0, : 1 : 0, and 1:0:0, and studied the theoretical and experimental 
expectations for the muon- neutrino flux at Earth, (j)^, and for the flavour ratios at Earth, 
R = <t>n/<f>e and S = cpr/cj)^. CPT violation (CPTV) has been implemented by adding an 
energy-independent term to the standard neutrino oscillation Hamiltonian. This introduces 
three new mixing angles, two new eigenvalues and three new phases, all of which have 
currently unknown values. We have varied the new mixing angles and eigenvalues within 
certain bounds, together with the parameters associated to pure standard oscillations. Our 
results indicate that, for the models 1:2:0 and : 1 : 0, it might possible to find large 
deviations for R, and S between the cases without and with CPTV, provided the 
CPTV eigenvalues lie within 10~ 29 — 10~ 27 GeV, or above. Moreover, if CPTV exists, 
there are certain values of R and S that can be accounted for by up to three production 
models. If no CPTV were observed, we could set limits on the CPTV eigenvalues of the 
same order. Detection prospects calculated using IceCube suggest that for the models 
1:2:0 and 0:1:0, the modifications due to CPTV are larger and more clearly separable 
from the standard-oscillations predictions. We conclude that IceCube is potentially able 
to detect CPTV but that, depending on the values of the CPTV parameters, there could 
be a mis-determination of the neutrino production model. 
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1. Introduction 

Experiments performed over the last thirty years have established that neutrinos can change 
flavour: careful measurements of solar fl|,|l, H]> atmospheric |||, reactor || and accelerator 
0, ||] neutrinos have established that there is a nonzero probability that a neutrino created 
with a certain flavour is detected with a different one after having propagated for some 
distance, and that this probability is a periodic function of the propagated distance, L, 
and the neutrino energy, E. The standard mechanism of neutrino flavour change requires 
neutrinos to be massive and results in a probability of flavour change that is oscillatory, 



with oscillation lengths that have a distinct 1/E dependence (see Section 2.1). 

So far, the experiments that have studied neutrino flavour transitions || have been 
designed to detect neutrinos with energies that range from a few MeV (solar neutrinos) 
to the TeV scale (atmospheric neutrinos). Notably, data from the Super-Kamiokande 
atmospheric neutrino experiment |l(| was used to find an energy dependence of the oscil- 
lation probability of E n , with n = —0.9 ± 0.4 at 90% confidence level, thus confirming the 
dominance of the mass-driven mechanism in this energy range, and relegating any other 
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potential mechanisms to sub dominance. It is possible, however, that one or more of such 
subdominant mechanisms become important at higher energies. 



In the present paper, we have explored a possible scenario where there is an additional 
oscillation mechanism present which results in an energy-independent contribution to neu- 
trino oscillations. This mechanism, though subdominant in the MeV-TeV range, might 
become dominant at higher energies, where the 1/E dependence of the standard oscillation 
term might render it comparatively unimportant: the higher the energy, the stronger the 
suppression of the standard oscillation term. The highest-energy flux of neutrinos available 
is the expected ultra- high-energy (UHE, with energies at the PeV scale and higher) flux 
from astrophysical sources -notably, active galaxies (see Section ||)- which are located at 
distances in the order of tens or hundreds of Mpc. 

We will deferr the detailed treatment of how the energy-independent contribution is 
introduced to Section || and focus now on the possible mechanisms that might motivate 
it. According to the CPT theorem, any Lorentz-invariant local quantum field theory must 
be built out of a CPT-conserving Lagrangian. However, the Standard Model (SM) is 
known to be valid at energies well below the Planck scale, mpi — 10 19 GeV and, at higher 



energies, motivated by theories beyond the SM [11, 12], CPT and Lorentz invariance might 
be broken. At accessible energies, the breaking of these symmetries can be described by 
an effective field theory that contains the SM. We have explored the possibility that CPT 
is not an exact symmetry, but rather that it is broken by the addition of a CPT-odd 
term to an otherwise CPT-even Lagrangian. The observation of the non-conservation of 
CPT would imply a fundamental revision of the usefulness of local quantum field theories as 
accurate descriptions of fundamental interactions. A possible realisation of a CPT-violating 
(CPTV) effective field theory is the Standard Model Extension [T3|] , which contains the 
SM, conserves SU (3) x SU (2) x U (1), and also considers potential Lorentz- and CPT- 
violating couplings in the gauge, lepton, quark and Yukawa sectors. It is worth noting that 
an alternative mechanism that also results in an energy-independent contribution to the 
oscillations is the nonuniversal coupling of the different neutrino flavours to an external 



torsion field 15]. 



The rest of the paper is organised as follows. In Section [||, we explain how the contri- 
butions from potential energy-independent new physics affect the flavour-transition proba- 
bility and introduce the detected fluxes of the different flavours of astrophysical neutrinos, 
and the ratios between them, as our observables. Section || discusses the deviations of 
the flavour ratios from their standard values in three different cases: when only can 
be detected, when and v e can be detected, and when the three flavours, u e , z^, u T , 
can be detected 1 . In Section we predict the feasibility of detecting a possible energy- 
independent contribution using the IceCube neutrino detector, and larger-volume versions 
of it. We summarise our results and conclude in Section [|. 



1 Some preliminary results on the three-flavour detection case were presented by the authors elsewhere 
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2. Theoretical framework 



2.1 Standard mass-driven flavour oscillations 

The standard mechanism that explains neutrino flavour transitions makes use of two differ- 
ent bases: the basis of neutrino mass eigenstates, which have well-defined masses, and the 
basis of neutrino interaction states -the flavour basis- which are the ones that take part in 
weak processes such as W decay. The two bases are connected through a unitary transfor- 
mation, so that we can write each one of the flavour states \v a ) as a linear combination of 
the mass eigenstates i.e., 

i 

where the coefficients [i7o] a j are components of the unitary mixing matrix that represents 
the transformation. Assuming the existence of three active neutrino families (a = e,/i,r), 
as indicated by LEP (TtJ, three mass eigenstates are required (i = 1,2,3 in the sum) and 
so Uo is a 3 x 3 matrix. The mass eigenstates |z/j) satisfy Schrodinger's equation and so 
propagate, sans a flavour-independent common phase e~ tEL , as 

\u { (L)> = e- iHL \ui) = e-^ L h) , (2.2) 



where we have assumed that m, <C E, so that p = W 1 E 2 — m? ~ E — m 2 /(2E), and 
that, because neutrinos are highly relativistic particles, t ~ L (in natural units). Thus, 
a neutrino created with a definite flavour a will, in general, become a superposition of 
states of different flavour as it propagates, so that a detector placed in its way can, with a 
certain probability, register it as having a different flavour from the one that it was created 
with. This becomes evident by writing the Hamiltonian in the basis of flavour eigenstates, 
a choice which will also allow us to introduce contributions from new physics later in a 
more straightforward manner. In this basis, if a neutrino is produced with flavour a, then, 
after having propagated for a distance L, its evolved state will be 

K(L)) = e - iH ™ L \u a ) , (2.3) 

where the oscillation Hamiltonian H m is the one corresponding to the standard, mass- 
driven, mechanism, and is written in the flavour basis. H m is related to the Hamiltonian 
in the mass basis -the "mass matrix"- through a similarity transformation that makes use 
of the unitary mixing matrix Uq: 

± diag (0, Arrini , Ami ) + , 
H m = U HUt = U &l 2g 21 ' ^-U\ . (2.4) 

Uq is the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix, which in the PDG 



parametrisation IS] can be written in terms of three mixing angles, #12, #13 and 623, and 
one CP-violation phase, 6cp, as 

/ C12C13 S12C13 si 3 e~ i5cp 



-S12C23 - Cl2S23Sl3e" 5cP C12C23 - Si 2 S 2 3Sl3e lScp S23C13 
\ S12S23 - C 12 C23Sl 3 e lScp -C12S23 - S 1 2C23Sl 3 e l5cp C23C13 
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with Cij = cos (Oij), Sij = sin (9ij). The standard flavour-oscillation probability P a p = 
\{vfj\v a (L)}\ 2 can hence be calculated, using Eq. (|2.3| ) for the evolved neutrino state and 
Eq. (|2.4|) for the standard Hamiltonian, to be 



P a , = ^ - 4 J> (4) sin 2 (^l) + 2 £ Im (j* ) sin (^l) , (2.5) 

i>j V / i>j \ J 

where Am 2 .- = m| — m|, with rrn the mass of the i-th. eigenstate, and 

J% = Mi P7oU [^o] aj [Efofo • (2-6) 



(For a detailed deduction, see, e.g., [IS].) It is straightforward to conclude from Eq. 



that flavour transitions occur because neutrinos are massive, particularly, because different 
mass eigenstates have different masses (clearly, if Am| = 0, no transitions occur), and 
because flavour states are not mass eigenstates. Note the 1/E dependence on the energy 
associated with this standard, mass-driven, oscillation mechanism. Note also that the full 
form of the PMNS matrix includes two extra Majorana CP-violation phases, a% and 0:2, 
which are identically zero if neutrinos are Dirac, so that the complete matrix is given by 

U x diag (V Ql/2 , e iQ2/2 , l) . (2.7) 

However, these phases do not affect the oscillations, so we have not included them in the 
definition of the standard mixing matrix Uq. 

Because, as was mentioned in Section [l], we will be considering UHE neutrinos of 
extragalactic origin, the flavour-transition probability in Eq. (^^) oscillates very rapidly, 
and we use instead the average probability, which is obtained by averaging the oscillatory 
terms in the expression, thus yielding 

(P a p) = ^ U °U 2 \lUoU 2 ■ (2.8) 

i 

When using the average probability, the information in the oscillation phase, including 
any potential CPTV energy-independent contribution, is lost. The modifications to the 
mixing angles due to CPTV, however, will still be present in the averaged version of the 
probability. 

Using the latest data from solar, atmospheric, reactor (KamLAND and CHOOZ) and 
accelerator (K2K and MINOS) experiments, the authors of |^0| performed a global three- 
generations fit and found the best-fit values of the standard oscillation parameters and 
their 3a intervals to be 

Am 2 .! = 7.65±o:[jo * 10~ 5 eV 2 , (Amf^ = 2.40±g;|| x 1(T 3 eV 2 (2.9) 

sin 2 (0 12 ) = 0.304±g;gg| , sin 2 (0 13 ) = 0M±° -° f , sin 2 (0 23 ) = 0.50l°;^ . (2.10) 

The are no experimental values for 5qp presently. We have assumed a normal mass hier- 
archy, so that Am| 2 = Am^ — Am^. 
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2.2 Adding an energy-independent Hamiltonian 

The lepton sector of the Standard Model Extension contains the Lorentz-violating contri- 



butions [13] 



Cuv => - (a L )^ L a -fL p + -i ( CL ) p L a -f D v Lp , (2.11) 



1 , , ^ „•«■, 
2 1 

where L a is the usual SM lepton doublet and a, j3 are flavour indices. The first term is 
CPT-odd, with the coefficients having dimensions of mass, while the second term is 
CPT-even, with the cl dimensionless. Because we are interested in an energy-independent 
CPTV contribution, we have kept only the first term in the Lagrangian. In the neutrino 
sector, then, CPT violation can be introduced through an effective, model-independent, 



vector coupling of the form [ 21 



£cptv= Trhtf-fJ • (2-12) 

The vector V^^v^ is CPT-odd and tff are real coefficients, so £cptv * s CPT-odd, i.e., 
CPT (££p TV ) = — ^CPTV" When the effective energy-independent Hamiltonian associated 
to £cptv ^ s added to the standard mass-driven neutrino oscillation Hamiltonian, it modifies 
the energy eigenvalues and, as a result, the mixing matrix is modified as well. 

Motivated by the vector coupling in Eq. ( 2.12j ), and in analogy to the standard oscil- 



lation scenario, we can introduce an energy-independent contribution in the form of the 
Hamiltonian (also in the flavour basis) 

H b = U b diag (0,621,631) Ul , (2.13) 

where bij = bi — bj. Following |ffif|, we write the mixing matrix in this case as 

U b = diag (l, e^.e* 8 ) U {{9 Hj } , 5 b ) . (2.14) 

The mixing angles associated with this Hamiltonian are #612, #613 5 $623) and S b fills the role 
of <5cp in the standard Hamiltonian. The two extra phases, (p2 and fa, appear because, 
once the flavour states and the mass eigenstates have been related through Eq. (p.l|), the 
former are completely defined, and the two extra phases cannot be rotated away. 

H b is dependent on eight parameters -two eigenvalues (621, 631), three mixing an- 
gles (#6i2, $613, $623) and three phases (6 b , fa, fa)- whose values are currently unknown. 
There are, however, experimental upper limits |EE on 621, obtained using solar and Super- 
Kamiokande data, and on 632, obtained using atmospheric and K2K data: 

621 < 1.6 x 1(T 21 GeV , 632 < 5.0 x 10~ 23 GeV . (2.15) 

The full Hamiltonian, including standard oscillations and the energy-independent con- 
tribution, is then 

H f = H m + H b . (2.16) 

In Section ||, we saw that H m has been experimentally demonstrated to be the dominant 
contribution to the oscillations in the low to medium energy (MeV-TeV) regime: there 
are no indications of new energy-independent physics at these energies and accordingly 
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the limits on bij shown in Eq. (|2.15| ) were placed. Because of the 1/E dependence of 
H m , however, it remains possible that, at higher energies, where the contribution of H m is 
suppressed, the effect of a hypothetical energy-independent term Hb becomes comparable 
to it or even dominant. Such energy requirement is expected to be fulfilled by the UHE 
astrophysical neutrino flux (see Section |IJ) . 

We would like to write the flavour transition probability corresponding to this Hamil- 
tonian in a form analogous to Eq. ( |2.8| ). In order to do this, we need to know what is the 
mixing matrix Uf that connects the flavour basis and the basis in which Hf is diagonal. Us- 
ing basic linear algebra, this is achieved simply by diagonalising Hf, finding its normalised 
eigenvectors, and building Uf by arranging them in column form. The components of the 
resulting matrix are in general complicated functions of the standard mixing parameters 
({Oij}, I Am^ j, 5cp) and of the parameters of Hb {{Obij}, {bij}, 5b, 4>2, 4>z)- By comparing 
the mixing matrix obtained by diagonalisation of Hf with a general PMNS matrix, given 
explicitly by Eq. (2.5) with mixing angles Oij and phase 5f, we are then able to calculate 
how the effective mixing angles @ij vary with the parameters of Hf, and 6cp- Succintly 
put, we have 

Uf = U f ({%} , {9 Mj } , {Am?,} , {btj} , S C P, 5 b , A2, fos) = U ({©,,} , 5 f ) . (2.17) 



Note that we have not used a perturbative expansion in the bij, as in [21], to calculate 
U f . This was done in order to allow for the possibility that the new physics effects become 
dominant at high energies, a possibility that would be negated if we had assumed that 
the effects are small from the start. As a result, the functional forms of the 0^ and 5f, 
while calculated in a straightforward manner, result in lengthy expressions and, due to 
their unilluminating character, we have chosen not to present them here. 

Thus defined, Uf is 14-parameter function. However, the standard mixing parameters 
Am^i, Am^, 9\2, 6*i3 and 623 have been fixed by neutrino oscillation experiments (see 
Eq. ( |2.9| )). Additionally, in order to simplify the analysis, we have set the phases 5b = 
02 = ^3 = 0. The standard CP-violating phase 5cp has been allowed to vary in the range 
[0, 2ir] in some of our plots, but otherwise we have set it to zero as well. As a further 
simplification, we have made the eigenvalues of H m proportional to those of Hb, at a fixed 
energy of E* = 1 PeV, that is, 

2 



Am., 

bi , = A ^ , (2.18) 

%3 2E* 



with A the proportionality constant. The upper bounds on the bij, Eq. ( p.15 ), are satisfied 



for A < 10 4 . Standard, purely mass-driven oscillations are recovered when A = 0. Thus 
we are left with only four free parameters to vary: A, #612, #613 and ^23 (and <5cp, where 
noted). 

In analogy to Eq. d2.8[) , the average flavour transition probability associated to the 
full Hamiltonian Hf is then 

(P«e)=Y,\lUf]J 2 \lUf]pi\ 2 . (2-19) 
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We will use this expression for the flavour-transition probability hereafter. It is worth 
noting that, if the CPTV contribution were introduced instead through a modified energy- 
momentum relation, only the oscillation phase would be affected, and this information 



would be lost when the average probability was used in place of the oscillatory one [ 22 



2.3 Astrophysical neutrino flavour ratios 

We have seen that, in order for a potential energy-independent contribution to the flavour 
transitions to be visible, we would need to use the expected UHE astrophysical neutrino 
flux. As mentioned in Section [l], the sources of this flux, e.g., active galaxies, are located 
at distances of tens to hundreds of Mpc, so that the average flavour transition probability, 



Eq. ( 2.19 ), can be used. 



If, at the sources, neutrinos of different flavours are produced in the ratios <§r e (j>% 
then, because of flavour transitions during propagation, the ratios at detection will be 



E ( p p«)4> ( 2 - 2 °) 



p=e,/Lt,T 

for a = e,fi, r. Note that the ratio in Eq. ( |2.20| ) is the proportion of v a to the sum of 
all flavours detected at Earth. We will later denote the actual neutrino fluxes, in units of 
GeV -1 cm~ 2 s _1 sr _1 , by <3? a . Evidently, the initial flavour ratios depend on the astro- 
physics at the source, which is currently not known with high certainty, while the detected 
ratios depend also on the oscillation mechanism and could be affected by the presence 
of an energy-independent contribution at high energies. Thus, the reconstruction of the 
initial neutrino fluxes from the detected ones is a difficult task [25, 24, 26, 27, 28, 29]. 



The effect of CPTV on the flavour ratios of high-energy astrophysical neutrinos has been 
explored elsewhere literature: in p0| , for instance, a two- neutrino approximation was em- 
ployed and it was assumed that H m and are diagonalised by the same mixing matrix, 



i.e., that Quj = while, in [23], neutrinos and antineutrinos were treated differently due 
to CPTV. Ref. [^lj used a formalism similar to the one we have used, but applied it to 
long-baseline terrestrial experiments and, due to the lower energies involved, introduced 
the CPTV effects as perturbations. The main difference between the existing literature 
on the effects of CPTV on the flavour fluxes of high-energy astrophysical neutrinos and 
the present work is that we have not treated the CPTV contribution as a perturbation, 
but, rather, we have allowed for the possibility that it becomes dominant at a high enough 
energy scale. 

The most commonly used assumption for the initial neutrino flux [^6| considers that 
the charged pions created in high-energy proton-proton and proton-photon collisions decay 
into neutrinos and muons, which in turn decay into neutrinos: 

7r + -)• fj+i/p ->■ e+VeV^v^ , 7r~ -> {TV M -» e~V e y f ^ v , . (2-21) 



Such process yields approximately ^ : : = 1 : 2 : (see [27] for a more detailed 
treatment), where we have not discriminated between neutrinos and antineutrinos, as is the 
case with current Cerenkov-based neutrino telescopes. In the standard oscillation scenario, 
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Production 


Initial flux 


Std. detected flux 


^>std 


gstd 


mechanism 




^std . ^std . ^std 






Pion decay 


1:2:0 


1:1:1 


1 


1 


Muon cooling 


0:1:0 


0.22 : 0.39 : 0.39 


1.77 


1 


Beta decay 


1:0:0 


0.57 : 0.215 : 0.215 


0.38 


1 



Table 1: Standard values (without energy-independent new physics contributions) of the detected 
flavour ratios <j) a {a — e, fi, r) and of the ratios-of-ratios R, S, for the three scenarios of initial 
flavour ratios considered in the text. The detected ratios were calculated using the average flavour- 
transition probability in Eq. (2.8) with the central values of the mixing angles: sin 2 (#12) = 0.304, 
sin 2 (0 13 ) = 0.01, sin 2 (6> 23 ) = 0.50. We have defined R std = <jf^ d /<j)f d and S std = 4>f d /^ d . 



i.e., in the absence of an energy-independent contribution, plugging this initial flux into 
Eq. (|2T20| ), and using the best-fit values of the mixing angles, Eq. (|2.9| ), results in equal 



detected fluxes of each flavour, i.e., 4>f '■ n '■ 4>t ~ 1 : 1 : 1- 

In a related production process [27, [H], 32 1, the muons produced by pion decay lose 
most of their energy before decaying, so that a pure-f M flux is generated at the source, i.e., 
0g : (j)^ : 0^ = : 1 : 0. In the standard oscillation scenario, these initial ratios result in 
the detected ratios (f)f d : <^ td : 0f d w 0.22 : 0.39 : 0.39. Alternatively, a pure-i/ e initial 
flux, corresponding to (jP e : <p° u : 0° = 1 : : 0, produced through beta decay has been 



considered, e.g., in [27]. In this scenario, high-energy nuclei emmitted by the source have 
sufficient energy for photodisintegration to occur, but not enough to reach the threshold 
for pion photoproduction. The neutrons created in the process generate v e through beta 
decay. For these initial ratios, the resulting detected ratios, in the standard oscillation 
scenario, are c/>| td : (f) s ^ d : 0^ td « 0.57 : 0.215 : 0.215. The results are summarised in Table 
|l[ In the following sections, we will consider the possibility of observing the hypothetical 
energy-independent contribution of Hf, assuming that the initial ratios correspond to one 
of these three production scenarios. 

Using the detected flavour ratios, we have defined the ratios of ratios 

R = h. , S = ^ . (2.22) 

Their standard values, R std and »S std , i.e., those calculated in the absence of H^, are shown 
in Table |l| for the three choices of initial flavour ratios. Note that S std = 1 for any choice 
of initial ratios because the value of 623 used was its best-fit value vr/4, which yields equal 
detected fluxes of and v T due to maximal mixing. Deviations from this value result in 
S 7^ 1 p3| . In the following sections, when we allow to contribute, we will calculate the 
extent to which the values of R and S deviate from their standard values. 

Since we are considering neutrinos that travel distances of tens or hundreds of Mpc, 
neutrino decay is a possibility. Flavour mixing and decays in astrophysical neutrinos have 



been explored before, e.g., in [34, 23] The current strongest direct limit on neutrino 
lifetime, r/m > 10 -4 s eV -1 , was obtained using solar neutrino data [p6| . More stringent, 
though indirect, limits can be obtained by considering neutrino radiative decays and using 
cosmological data [37]: r > few x 10 19 s or t > 5 x 10 20 s, depending on the mass hierarchy 
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and the absolute mass scale. Assuming that the heaviest mass eigenstates decay into the 
lightest one plus an undetectable light or massless particle (e.g., a sterile neutrino), then, 
following p4l], the ratio of flavour a at Earth will be 



0=e,/j,,T i (3=e,/i,ri(stable) 

(2.23) 

where r, is the lifetime of the i-th mass eigenstate in the laboratory frame. As explained 



in [24], this expression corresponds to the case where the decay has been completed when 



the neutrinos arrive at Earth. In a normal hierarchy, v\ is the only stable state and so 



^dec,norm 



WoU 2 E ( 2 - 24 ) 



while in an inverted hierarchy 1/3 is the stable state and 

<e c ' inv = l[^U 2 £ ^|[£/oW 2 ■ (2-25) 

/3=e,/i,T 

Reference |24|] provides expressions for the flavour ratios when the decay product accom- 
panying the lightest eigenstate can also be detected, but since the focus of our analysis 
is the modification of the oscillations through terms of the form Eq. ( [2.12| ), we have cho- 
sen to include only the simplest case of neutrino decay, described by the two preceding 
expressions. 

The reader should bear in mind that there are other mechanisms involving new physics 
that might also affect the astrophysical flavour fluxes at Earth, such, non-standard inter- 
actions [38], coupling of neutrinos to dark energy [39], deviations from the unitarity of the 
PMNS mixing matrix [40], to name just a few. 



3. Effect of an energy-indepedent contribution on the astrophysical flavour 
ratios 

3.1 Detection of 

In Figure [l] we present a plot of the muon- neutrino flavour ratio, (j)^ as a function of 
A. The coloured bands correspond to different neutrino production scenarios, namely: 
: (()® : (f>® = 1 : 2 : (blue), 0:1:0 (purple) and 1:0:0 (brown), which have 
been generated by varying the three CPTV angles (0&12, #613, #623) within [0,7r], and 
the three standard mixing angles (612, #13, #23) within their 3a bounds, with <5cp = 0. 
For comparison, we have included the hatched bands which represent the pure standard 
oscillation case, that is, without CPTV. These have been generated by setting A = and 
varying the three standard mixing angles within their 3a bounds, again with <5cp = 0. As 
expected, the standard-oscillation bands are contained within the corresponding CPTV 
region. 

When CPTV is allowed, we observe large deviations of (p^ with respect to the pure 
standard-oscillation bands, especially for the scenarios 4>® : (f>® : <^>° = 1 : 2 : and 0:1:0, 
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Figure 1: Allowed regions of values of the detected muon- neutrino flavour ratio, </> M , as a function 
of A for different neutrino production models. The three standard mixing angles (612, #13, #23) were 
varied within 3c bounds and the three CPTV angles (#612, &bi3, $623) were varied within [0,7r], 
while Sep = 0. The hatched regions are the allowed regions of </> M when only standard oscillations 
are allowed, and allowing the 9ij to vary within their 3er bounds. 



and less so for the scenario 1:0:0. Starting from A ~ 0.1 (632 - 1.2 x 10~ 28 GeV), the 
CPTV bands start growing with A, as was expected, since A measures the strength of the 
CPT violation. Thus, the influence of the CPTV contribution to the oscillations grows 
and, as a consequence, the accessible region also grows. This is due to the wide range of 
values that the CPTV mixing angles can take, in comparison with the standard ones. Past 
A = 1, the CPTV regions reach a plateau, owing to the fact that the CPTV term becomes 
dominant over the standard-oscillation term in the Hamiltonian. 

An interesting feature is the overlap between the standard-oscillation band for the sce- 
nario 1:2:0 and the CPTV region for scenario 0:1:0. A similar overlap occurs between 
the scenarios 0:1:0 (without CPTV) and 1:0:0 (with CPTV). As a consequence of 
these overlaps, if CPTV exists for certain values of the parameters, a measurement of 
will be insufficient to distinguish what the neutrino production model is. For instance, if a 
value of 4>^ ~ 0.4 were measured, and A > 0.2, we would not be able to assert whether the 
initial fluxes were : 1 : or 1 : : 0. Analogously, since the standard-oscillation bands are 
contained within the CPTV regions, for these cases we will be unable to conclude, from 
the measurement of </>«, whether or not CPTV effects are present. 

Although we have not presented it here, we have tested that when varying 5&P within 
[0, 27r], the regions change only in a few percent. Therefore, the features observed in Fig. |j] 
are largely independent of the value of Sep- 
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Figure 2: Allowed regions of values of i? = ^/0 e as a function of A, for three scenarios of initial 
flavour ratios: <jP e : <jp : <\P T — 1 : 2 : (blue), 0:1:0 (purple) and 1:0:0 (brown). For each one 
of them, at each value of A, the three CPTV mixing angles, #bi2, Qb\z and Qbiz, were independently 
varied within [0, 7r], and the three standard mixing angles, #12, #13 and #23, were varied within their 
3<t bounds (Eq. ( |2.9D ): the lowest and highest value of R obtained in this way define, respectively, 
the lower and upper bounds of the corresponding region, for this particular value of A. The CP- 
violating phase 5qp = 0. The upper horizontal axis is 632 = A Am| 2 / (2E*), with E* = 1 PeV; 
to find the corresponding values of 621, notice that &21/&32 = Am^/Aml ~ 1/30. The hatched 
horizontal bands are the allowed regions of R assuming only standard oscillations and allowing the 
standard mixing angles 9ij to vary within their 3cr bounds. 

3.2 Detection of two flavours: and v e 

In this section, we consider an scenario where only muon- and electron-neutrinos are de- 
tected in an available neutrino telescope. Assuming that flavour identification is possible, 
we can define the ratio of observed flavour ratios 

R = <j)J<j) e , (3.1) 

which, in the presence of CPTV, depends on A and on the three CPTV mixing angles 
through the definition of the <p a in terms of the flavour-transition probabilities. Figure 
|2] shows the regions of values of R as a function of A by varying the three CPTV angles 
within [0, 7r]. 

The coloured regions in Fig. || grow with A, for the same reason as they did in Fig. [l]. 
Under the assumption of a : 1 : production model, R can attain very large values, 
between 10 6 and 10 7 , as a result of very low electron-neutrino fluxes. On the other hand, 
the CPTV region associated to the 1:2:0 production model reaches a maximum of R = 2 
after A ~ 5, while the one associated to 1 : : reaches a maximum of R = 1 after A ~ 1. 
Therefore, a measurement of R 4 could imply that the production model is : 1 : 0, 
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Table 2: Conclusions that can be obtained depending on the measured value of R, according 
to Fig. H [i?i2o] > [-Roio] an d [-^ioo] represent, respectively, the standard-oscillation bands corre- 
sponding to the <jP e : <jp : (f) T — 1 : 2 : 0, : 1 : and 1:0:0 production models. When we refer 
to standard oscillations, we mean either the inexistence of CPTV (A = 0) or the existence of too 
small a CPTV (A < 1). 



and that CPTV effects are present, but will not be enough to set strong bounds on A. The 
minimum value in both the 0:1:0 and 1:2:0 regions is located around R = 0.8 ~ 1, 
while for the 1:0:0 flux, it is zero for most of the range of A. 

If a value of R < 4 is found, the ability to single out a production model depends on 
the exact value of R that is measured. A single production model can be distinguished 
univocally for some measured values of R (e.g., for R < 0.84, the model is 1 : : 0), while 
for others, two (when 1 < R < 4, depending on A) or three models (when R ~ 1) can 
account for the same measured value. In the same way, for some values of R, CPTV and 
standard oscillations cannot be distinguished. The conclusions, from Fig. ^, that can be 
obtained from the measurement of different values of R are shown in Table [|. 

3.3 Detection of three flavours: u^, v e and v T 

Given the approach of this work, we have considered as a natural step to extend our analysis 
by taking into account the possibility of tau-neutrino detection. For this purpose, we have 
defined the ratio 

S = (j) T /<f)^ , (3.2) 

and studied the effects of the new physics in the R vs. S plane. 

Before we describe the results of this section, it will be useful to turn our attention to 
the values of R and S associated to the different production models when only standard 
oscillations are allowed, along with the correspondong standard values of the detected 
fluxes 4> a . These values are shown in Table [l]. The CPTV regions in Figures [3|-|5] have 
been generated by setting A = 100 (Figs. ||] and f|) and A = 1, 10, 100 (Fig. |5|), fixing the 
standard mixing parameters at their best-fit values (see Eq. (|2.9[) ), and varying the CPTV 
mixing angles O^j within the range [0,7r], while the standard-oscillation regions have been 
generated by setting A = and varying the standard mixing angles within their 3<r bounds. 
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Figure 3: Regions of R and S" accessible with CPTV by assuming different neutrino production 
scenarios. Different colours correspond to different initial ratios: <fP e : 0" : </>° = f : 2 : (blue), : f : 
(purple), f : : (brown). Darker regions are generated with standard (CPT-conserving) neutrino 
flavour oscillations, by allowing the standard mixing angles to vary within their 3cr experimental 
bounds. Lighter regions correspond to the case when we include a dominant CPTV contribution 
with A = fOO, and allow the CPTV angles Ouj to vary within [0,7r]. All phases are set to zero. 

With the exception of Figure ^, where <5cp has been allowed to vary, we have set all phases 
equal to zero. 

In Fig. |3| we display in three R vs. S panels the allowed regions of values that correspond 
to pure standard oscillations at 3a, in dark tones, and the corresponding regions allowed by 
the mixed solution composed of the standard oscillation plus CPTV effects, with A = 100, 
in lighter tones. At this value of A, the CPTV contributions are the dominant ones in 
the Hamiltonian, i.e., Hf ~ H^. Each plot corresponds to a different neutrino production 
scenario: (a) to 1 : 2 : 0, (b) to : 1 : 0, and (c) to 1 : : 0. 

From these plots we can extract two observations: one is the potentially dramatic 
deviation of the allowed values of the pair (R, S) when CPTV is turned on, and the other 
is the presence of points that are common to all scenarios. The latter implies that, if 
there is CPTV, there could exist pairs (R,S), such as (1,0), that can be generated by any 
of the three production models, each with a different set of values for the CPTV mixing 
angles. There are also (R, S) pairs which could be generated by one production model 
with standard oscillations or with a different model with CPTV, e.g., those lying around 
(1,0.95). It is convenient to remark that this figure and Fig. |2| are consistent with each 
other, which can be shown by projecting the CPTV regions of Fig. ^ onto the horizontal 
axis and checking that the limits on R agree with those on Fig. |2|. We have marked a few 
notable points in each plot: for the three production models, the points labelled with A 
correspond to the best-fit values of the standard-oscillation mixing parameters in Eq. (|2.9|). 

In Fig. |], the allowed R — S regions corresponding to the three neutrino production 
models are shown together: 1 : 2 : in blue, : 1 : in purple, and 1 : : in brown, 
where, as before, the darker tones correspond to pure standard oscillations, and the lighter 
tones, to a dominant CPTV with A = 100. Following the argument in Section |2.3| , we have 
also included the regions of (R, S) pairs allowed by neutrino decay into invisible products, 
considering both the cases of a normal mass hierarchy (decay into the v\ eigenstate), in 
red, and of an inverted hierarchy (decay into z^), in green. The decay regions were also 
generated by varying the standard-oscillation mixing angles within their 3a bounds. 



-13- 



1.4 



1.2 



1.0 



0.6 



0.4 



0.2 



0.0 



R = <l>^l<l>e S = 4> T I4>^ 




■ | 1:2:0 std. osc. (A = 0) 
I | 1 :2:0withA = 100 
0:1:0 std. osc. (A = 0) 
0: 1 :0withA= 100 



| | 1:0:0 std. osc. (A = 0) 
H l:0:0withA= 100 

decay to V\ 
^ decay to V3 



10 



15 



Figure 4: Regions of R and S accessible by assuming different neutrino production scenarios: 
(jP e : 0° : 0° = 1 : 2 : (blue), : f : (purple) and 1:0:0 (brown). Darker shades of blue, purple 
and brown correspond to standard, CPT-conserving (i.e., A = 0), flavour transitions, while lighter 
shades correspond to flavour transitions dominated by CPT violation (A = 100). The former were 
generated by varying the standard mixing angles dij within their current 3cr experimental bounds; 
the latter, by fixing the 9^ to their best-fit values and varying the new mixing angles Ouj within 
and 7T. All of the phases were set to zero: i5cp = <5& = <fib2 = <t>b3 — 0. Also shown are the 
regions accessible through decay of the neutrinos into invisible products, when v\ is the lightest 
mass eigenstate (red) and when is the lightest one. 



While the regions corresponding to 1 : 2 : do not overlap those corresponding to 
neutrino decay, this is not the case for the 0:1:0 scenario, for which there is a clear 
overlap between the CPTV region and the one from neutrino decay into v%. The region 
corresponding to decay into v\ is completely contained within the 1:0:0 CPTV region. 
Projecting onto the horizontal axis, we see that, apart from the same overlaps between 
the allowed R intervals that existed in the two-neutrino CclSG, ELS observed in Fig. new 
overlaps appear due to the inclusion of the neutrino-decay regions. Particularly, whereas 
prior to the inclusion of decays, a value of R 3> 4 was a clear indication of the existence of 
CPTV and of a : 1 : production model, now Fig. |I| shows that, after R = 7, the decay 
into z^3, assuming an inverted mass hierarchy, can also be accountable for a high value of 
R. On the other hand, values of 0.25 < R < 0.75 can be reached by the 1:0:0 production 
model, with or without CPTV, as well as with neutrino decay into u%, assuming a normal 
hierarchy. A normal mass hierarchy is able to yield values of S as low as ~ 0.075, leaving 
only the small window between zero and this value as unique feature of CPTV. Taking into 
account experimental uncertainty, it is likely that this window is actually non-existent. 
Only if the measured S > 1.1, and in the absence of decays, would it be possible to identify 
a single production model, 1 : : 0, as the one responsible. If decays are allowed, however, 
they can also account, irrespectively of the mass hierarchy, for S > 1.1, and the ability to 
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Figure 5: Regions of values of R and S accessible when varying the parameter A between (no 
CPT breaking) and 100 (dominant CPTV term), for different neutrino production models. In blue: 
(jP e : 0° : (j)® = 1 : 2 : 0; in purple, : 1 : 0; in brown, 1:0:0. The region corresponding to pure 
standard oscillations (A = 0) was generated by varying the standard mixing angles within their 
3(J experimental bounds, Eq. ( |2.9| ). To generate the regions corresponding to A = 1, 10, 100, the 
standard mixing angles were fixed to their best-fit values, and the three CPTV mixing angles 9uj 
were varied independently within [0,7r]. The CP-violating phase <5cp = for all of the regions. 



single out a production model is lost. Thus, the signatures of CPTV become less unique 
in the presence of decays. A more complete analysis of neutrino decays [23], exploring also 
the possibility of incomplete decays and decay into visible products, further reduces our 
ability to uniquely identify the presence of CPTV. 

In the absence of neutrino decays, a measurement of R > 4.1 and S < 1.1 could 
indicate the presence of dominant CPTV and a : 1 : production model. Also, for 
R < 0.9, the production model is 1 : : 0; if S < 0.45, there is dominant CPTV and, for 
other values of S, the existence of CPTV will depend on the value of R measured. Close 
to R = 1, and for S < 0.9, any of the three production models with dominant CPTV can 
explain the measured (R, S) pair, while for 0.9 < S < 1, the pair could be generated either 
by the : 1 : or 1 : : models with CPTV, or by the 1:2:0 model with standard 
oscillations. For 1.1 < R < 4 and S < 0.6, there is a large region of overlap between the 
1:2:0 and 0:1:0 models with dominant CPTV. 

Fig. U shows the effect of the variation of the parameter A between (no CPTV, i.e., 
Hf = H m ) and 100 (dominant CPTV term, i.e., Hf ~ Hb) on the allowed R — S regions, 
for the different neutrino production scenarios. For the three of them, we can observe 
significant deviations from the predictions of the standard-oscillation case, even for low 
values of A. In this sense, it is interesting to point out that in the case of an experimental 
non-detection of CPTV in the neutrino flavour ratios, R and S can be used to set limits 
on the related parameters. In fact, when A = 1, and for a neutrino energy of 1 PeV, we 
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Figure 6: Variation of R and S with $13, when the CP-violation phase 5qp is allowed to vary 
between and 2tv. The upper limit for #13 is given by the current bound sin 2 (#13) < 0.056 (3<r) 
and the mixing parameters #12, #23, Am^ and Am| 2 have been set to their current best-fit values; 
standard flavour oscillations have been assumed throughout (i.e., A — 0). Three different scenarios 
of initial flavour ratios have been considered: <fP e : 0° : 0° = 1 : 2 : (in blue), : 1 : (in purple) 
and 1 : : (in brown). To allow for comparison, the arrows point to the curves on which <5cp = 0. 



can attain limits for the CPTV eigenvalues bij in the order of 10 and 10 GeV, for 
621 and 6237 respectively. It is also important to mention that these results can be easily 
rescaled to any energy, just by doing bijX(PeV/E). This would mean a very significant 
improvement over the current bounds of 10~ 23 — 10~ 21 GeV for 621 and 632, respectively 




Finally, it is worth exploring the effect of the choice of value of <5cp on R and S. From 
Fig. |(| we see that the effect on R of a non-zero value of 5cp in the standard-oscillation 
regions is more prominent for a choice of initial ratios of : 1 : 0, and less so for 1:2:0 
and 1:0:0. For S, the effect of a non-zero phase is greater for 1:0:0, and less for 0:1:0 
and 1:2:0. Note that using a non-zero value of Sqp can lead to a variation in R of up to 
42% (in the 0:1:0 case) and in S of up to 40% (in the 1:0:0 case) with respect to the 
standard case of 6*13 = 5qp = 0, when the largest 3a allowed value of sin 2 (#13) = 0.056 is 
assumed. 



4. Detection prospects at IceCube considering CPTV 
4.1 Experimental setup 



The IceCube neutrino telescope [41, 42 1, located at the South Pole, can detect both muons 
and showers initiated by incoming high-energy astrophysical neutrinos. Muon- neutrinos are 
detected through the high-energy muon produced in charged-current (CC) deep inelastic 
neutrino-nucleon scattering: the Cerenkov light emitted by the fast-moving muon in ice 
is detected by the photomultipliers buried in ice and used to reconstruct the muon track. 
Electron-neutrinos can generate electromagnetic and hadronic showers in CC interactions. 
Tau-neutrino CC events have a distinct topology: they either show two hadronic showers 
joined by a tau track (a double bang) or a tau track ending with the decay of the tau in 
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a hadronic shower (lollipop). Additionally, every flavour of neutrino is able to generate 
hadronic showers in neutral-current (NC) interactions. Flavour identification is expected 



to be difficult at IceCube [24| and the number of lollipops expected is about one every two 
years (for a neutrino flux of 10 -7 GeV -1 s _1 cm" 2 sr -1 [24Q, so the useful observables turn 
out to be the total number of NC plus CC showers, iV s h = -/V^ c + and the number 
of muon tracks, N v . From this, we can construct the closest experimental analogue of 
the variable R = 4>^/4> e as -R ex p = N v /N^. Due to the very low number of tau-neutrinos 
expected at IceCube, there is no practical experimental analogue of S = (pr/cp^. 

Following our analysis of CPTV in the previous sections, where we set the scale of 
CPTV atE* = l PeV, we have adopted the energy range Ef m = 10 6 < E v /GeV < £™ ax = 
10 12 for our predictions. Within this energy range, the Earth is opaque to neutrinos due 



to the increased number of NC interactions which degrade their energy [45|, so we have 
calculated only the number of downgoing events. Bear in mind, however, that due to the 
tighter background filtering that is required for the observation of downgoing neutrinos, 
which we have not taken into account in the following calculations, our estimates may be 
optimistic. Conveniently, for most of this range the atmospheric i/„ background flux will 
be well below the fluxes from the neutrino production models that we have probed (see 
Figure 0). Using the expressions of p3| , we can estimate the number of CC and NC events 



at IceCube, for a given astrophysical all-flavour diffuse neutrino flux, denoted here by 3> 



fall. 
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where E' v is the energy of the secondary neutrino in the NC interaction, T is the exposure 
time, tit = 5.1557 x 10 23 cm -3 is the number density of targets (nucleons) in ice, V e g is the 
effective detector volume (1 km 3 for IceCube) and f2 = 5.736 sr is the detector's opening 
angle (up to 85°). The total cross sections for neutrino and anti-neutrino deep inelastic 



scattering, (Jqq, c^C' ^nc an( ^ °NC' nave been extracted from [44]. The differential cross 
sections used are written as a function of both the primary and the secondary neutrino 
energies, which in NC interactions are related through E' v = (1 — (?/nc)) E v , with (j/nc) 
the NC inelasticity parameter, extracted from ||45|| . If the interval of interest (10 6 — 10 12 
GeV) is partitioned into small enough subintervals, then within each, the NC cross section 
can be approximated by ctnc = -AJSj? with A, B constants for each subinterval. We can 
then write 



vN _ a TpB _ a 



B 



(l-^NC^)))" 1 ^ , (4.3) 



so that 



<la ^ (E v , E' u ) = AB [1 — (y NC (E u ))]- B (Elf' 1 . (4.4) 



dE v 
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Using this expression in Eq. (4.2) and performing the E' v integral, we obtain the simplified 
form 

£max 

~ \ Tn T V eS Q / " dE v W* {E v ) x 

sh 

°NC Ok ) ~ °NC ~ ^sh ) , °NC Ov ) ~ ^NC Of ~ ^sh 



(4-5) 



with S and B taking the appropriate values in each subinterval. 

The number of CC showers generated by electron- and tau-neutrinos are, respectively, 
Ng° e = &JVg° and Ng c T = faN™, where 

€[0,1] (4.6) 



Pe + 9ti + \ 

are the normalised flavour ratios. The number of CC showers is given by 

Ng C = Ng c e + Ng c T = (4> e + <^ T ) (4.7) 
and the total number of CC plus NC showers is therefore 

N sh = Ng c + Ng c = (l - 4>,) N™ + N™ , (4.8) 

where we have used the fact that ^2p =e ^ T (ftp = 1- In a similar way, the number of 
downgoing muon-neutrinos is given by 

Nu, = • (4.9) 



Using Eqs. Q, (gj), (|J) and (gj), we find 



^cx P = = t ^ , (4.10) 



and we see that R exp depends on only one of the normalised flavour ratios, (f>u, and that it 
is independent of the exposure time and the effective detector size. Increasing T and V e g, 
however, results in a larger event yield and consequently in lower statistical uncertainty. 
Considering iV£\r and N~*9 as independent variables with Poissonian errors, i.e., \/ 



and \l respectively, we find the error on R exp to be 



As expected, <x# cxp a (TV e s) 1 , so that, for a given neutrino flux, the statistical error on 
-R exp decreases with the time of exposure and the size of the detector. 
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Figure 7: High-energy astro-physical neutrino flux models as function of neutrino energy. The 
spectral index a = 2.7 for the Becker-Biermann (BB) and 2.3, 2.6 for the Koers-Tinyakov (KT) 
fluxes with and without source evolution, respectively. For the BB flux, we have set r„/TcR = 3 
and -z™R X = 0-03 (see p7| ). The atmospheric muon- neutrino flux, modelled according to ^8j, which 
is considered as a background to the astrophysical neutrino signal, lies below the predictions of 
these models in the plotted energy range. 



Flux 


Vail 


N cc 

fall 




Waxman-Bahcall Q 

Becker-Biermann a = 2.7 [E7[] 

Koers-Tinyakov no source evolution a = 2.6 [48 

Koers-Tinyakov strong source evolution a = 2.3 [48 


1781.64 
9248.58 
9013.86 
16495.92 


33.12 
130.68 
354.96 
1866.42 


53.79 
70.77 
25.39 
8.84 



Table 3: Expected number of neutral-current and charged-current events (summed over all 
flavours) at IceCube, and N^, respectively, in the energy range 10 6 < E v /GcV < 10 12 , 

for different choices of the incoming astrophysical neutrino flux. The exposure time used was 
T = 15 yr, and the effective detector volume, V e g = 1 km 3 . Only downgoing events are consid- 
ered. The Waxman-Bahcall flux assumes an E~ 2 spectrum, while for the Becker-Biermann and the 
Koers-Tinyakov fluxes we have used a power law of the form E~ a , with a specified for each model. 
The different values of a have been selected so that the upper bound on the diffuse astrophysical 
muon-neutrino flux set by AMANDA-II is satisfied. Details of the neutrino production models 
can be found in the indicated references. 



4.2 Astrophysical neutrino flux models 

We will present our estimations of iV~9, and JV^/iV^, in the context of four 

models of astrophysical neutrino fluxes. The Waxman-Bahcall [46| model makes use of 
the observation of ultra- high-energy (> 10 19 eV) cosmic rays to set an upper limit on 



the neutrino flux. The limit depends on the redshift evolution of the neutrino sources, 
which could be active galactic nuclei (AGN) or gamma-ray bursts. We have adopted, 
conservatively, 

$wb = 1Q -8 (^/GeV)" 2 GeV" 1 cm" 2 s" 1 sr" 1 . (4.12) 

The second model, by Becker-Biermann [47], describes the production of neutrinos in the 
relativistic jets of FR-I galaxies (low-luminosity radio galaxies with extended radio jets) 
through the decay of pions produced in the interaction of shock-accelerated protons with 
the surrounding photon field. The sources are assumed to evolve with redshift according 
to certain luminosity functions. The flux is given by 

$™ (E u ) ~ 5.4 x 10" 3 (£„/GeV)~ 2 - 7 GeV" 1 cm" 2 s" 1 sr" 1 . (4.13) 



The third and fourth models, by Koers-Tinyakov [48], predict an astrophysical neutrino 
flux based on the assumption that the neutrinos originate predominantly at AGN and that 
these behave like Centaurus A, the nearest active galaxy. The difference between these 
last two models lies in the assumption about the redshift evolution of the source: in one 
of them, sources are assumed not to evolve with redshift, while in the other one, they are 
assumed to have a strong redshift evolution, following the star-formation rate of ~ (1 + z) 3 . 
These two fluxes are given, respectively, by 

$KT, no cvol. „ 3 5 x 1Q -10 (^/QeV) " L6 X 

x mm(l,E u /E U)hT ) GeV" 1 cm" 2 s" 1 sr" 1 (4.14) 
$ kt, cvol. _ 46 x 1Q -i2 (^/GeV)- 1 - 3 x 

xmm(l,E u /E u>hT ) GeV" 1 cm" 2 s" 1 sr" 1 , (4.15) 

with E v br = 4 x 10 6 GeV the break energy. Both the Becker-Biermann and Koers-Tinyakov 
models account for the change in cosmic-ray energy due to the adiabatic cosmological ex- 
pansion, but only the latter take into account also energy losses due to pion photopro- 
duction and electron-positron pair production in the interaction with the CMB photons. 
Currently, the most stringent upper bound on the astrophysical muon-neutrino flux is the 
one obtained by the AMANDA-II experiment [49], which restricts the integrated muon- 
neutrino flux to be lower than 7.4 x 10" 8 GeV cm" 2 s" 1 sr" 1 , within the interval 16 
TeV - 2.5 PeV. We have checked that the four fluxes used in our analysis satisfy this 
upper bound in the standard-oscillation case, i.e., when the detected flavour ratios are 
4> e : 4>u : <j) T = 1 : 1 : 1. The spectral indices of the power laws for the Koers-Tinyakov 
fluxes, —1.6 and —1.3, have been chosen so that the integrated fluxes between 16 TeV and 
2.5 PeV yield exactly the upper bound value set by AMANDA. A plot of the four different 
fluxes is presented in Fig. [?]. 

We have evaluated Eqs. ( |4.1| ) and (4.5) numerically in the range 10 6 < E^/GeV < 10 12 
for the four different fluxes, assuming T = 15 yr and V e Q = 1 km 3 (the IceCube effective 
volume). The results are presented in Table ||. The Waxman-Bahcall model yields the 
lowest number of CC and NC events, while the Koers-Tinyakov model with strong source 
evolution yields the highest number, more than two orders of magnitude over Waxman- 
Bahcall. 
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Figure 8: i? C xp = N^^/N^ vs. A for two neutrino flux models: the Waxman-Bahcall flux (in blue) 
and the Koers-Tinyakov flux with spectral index a = 2.3 and strong source evolution (in green). 
Three initial flavour ratios have been considered: </>g : <j>°^ : <fP T = 1 : 2 : (plots (a) and (d)), 0:1:0 
(plots (b) and (e)) and 1:0:0 (plots (c) and (f)). The coloured areas, bounded by dashed coloured 
lines, have been calculated plugging the values of M from Fig. [l] into Eq. ( 4.10| ) and using the event 
count data in Table ||. The dot-dashed coloured lines mark the la uncertainty on i? cxp for each 
assumption of initial ratios, calculated with Eq. (4.11). For comparison, the standard value of R oxp 
for the different initial ratios are shown in hatched bands of the corresponding colours; these values 
are calculated assuming no CPTV effects and using the current 3er bounds on the standard mixing 
angles. In the upper row, an effective detector area V^g = 1 km 3 (IceCube size) was assumed; for 
the bottom one, V e g = 5 km 3 (IceCube x 5). 



4.3 Results 

In this section, we will display the regions of ii xp vs. A, allowing for CPTV, for the 
Waxman-Bahcall model and the Koers-Tinyakov model with strong source evolution, which 
correspond, respectively, to the cases with the lowest and highest event yields. In order 
to generate these regions we have used the values in Table || to calculate R exp in the 
presence of CPTV, Eq. (4.10), using for (j)^ the range of values allowed by the variation 



of the standard mixing angles 9™ within their current 3<r experimental bounds and of 
the new mixing angles duj in the range [0, it] . This range of values of (p^ vs. A can be 
obtained from the plot of 4>n vs. A shown in Fig. |l|, after applying the transformation 
in Eq. (f4.6|). The resulting regions as functions of A are shown colour- filled in Figure || 
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for the Waxman-BahcaU (blue) and Koers-Tinyakov (green) models. We have explored 
two different detector effective volumes, 1 km 3 (IceCube-sized) and 5 km 3 , and the three 
different choices for the initial fluxes <$* e : ^ : (f)®. = 1 : 2 : 0, : 1 : and 1:0:0, that 
we introduced in the previous sections. Furthermore, this figure includes boundaries of la 
statistical uncertainty on R ex p that were obtained by adding (subtracting) la to (from) 
the upper (lower) boundaries of the coloured regions. The values of a were calculated by 
plugging into Eq. fl4.11| ) the corresponding values of cp^ that occur on the borderlines of 
the coloured regions. 

We note that the shapes of the coloured regions in Fig. |8| are similar to those in Fig. [I]. 
This can be understood if we note that R cxp , Eq. (|4.10| ), is proportional to In Fig. [j] 
we saw that there exists overlap among the regions related to the different hypotheses of 
the initial flux. In fact, in Fig. ^, there is overlap between the three production models, 
with the region corresponding to <fr® : 0^ : 0° = 1 : 2 : almost entirely enclosed within 
the region for 0:1:0. This last fact can be explained on account of the reduction in 
size of the region for 1:2:0, caused by the normalisation of (f>n applied for this case, i.e., 
(j)^ = 0^/3. As we can anticipate, due to the higher event yield, the statistical uncertainty 
on .Rexp associated to the Koers-Tinyakov flux is lower than the one associated to the 
Waxman-Bahcall flux, and is reduced when the larger detector volume is used. The size of 
this uncertainty is also proportional to the value of (see Eq. ( 4.11D ). As a consequence, 
the size of the la regions is larger for an initial flux of 1 : 2 : 0, intermediate for 0:1:0, 
and smallest for 1:0:0. 

When a detector volume of 1 km 3 is considered, there is a clear overlap among regions 
corresponding to different assumptions of the neutrino flux model when the production 
scenario is 1 : : 0. This observation is reinforced when we consider also the regions 
spanned by the statistical uncertainty. In comparison, for the 1:2:0 and 0:1:0 
scenarios, the regions corresponding to the two flux models do not overlap. When the 
5 km 3 detector is assumed, the regions associated to 1 : 2 : and 0:1:0 are further 
separated at the la level. However, there is still an overlap between the two flux models 
in the 1:0:0 scenario. Here we do not show the results for the Becker-Biermann model 
and the Koers-Tinyakov model with no source evolution, since their results are embodied 
in what we have already presented. For instance, the mean value of R exp and cx,R exp for 
Becker-Biermann is similar to the one for Waxman-Bahcall. Similar estimates can be easily 
obtained by plugging the values of Table ||| into Eqs. ( |4.10|) and ( [4.11 ). 



5. Summary and conclusions 

Motivated by the CPT-violating (CPTV) neutrino coupling considered in the Standard 
Model Extension, we have added a CPTV, energy-independent, contribution to the neu- 
trino oscillation Hamiltonian and explored its effects on the flavour ratios of the high-energy 
(1 PeV and higher) astrophysical neutrino flux predicted to come from active galactic nu- 
clei. We have parametrised the strength of the CPTV contribution by the parameter A, 
defined as the quotient between the eigenvalues of the CPTV Hamiltonian, 621 and 632, 
and those of the standard-oscillation one, Aml 1 /(2E*) and Am§ 2 /(2£*), with E* = 1 
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PeV, and allowed A to vary between 10~ 2 and 100, corresponding to standard-oscillation 
dominance and CPTV dominance, respectively. We have used three different neutrino pro- 
duction scenarios for the flavour ratios at the astrophysical sources: production by pion 
decay, which results in (fP e : 0° : 0° = 1 : 2 : 0; muon cooling, which results in : 1 : 0; 
and neutron decay, resulting in 1 : : 0, and explored the effect of a potential CPTV on 
the neutrino flavour ratios at Earth, <p a {a = e,fi,r), and on the ratios between them. 
With this objetive, we have studied the behaviour of 0^, R = (f>^/(f) e , and S = 
by letting the standard-oscillation mixing parameters vary within their current 3a experi- 
mental bounds, and varying the unknown CPTV parameters as broadly as possible, while 
keeping 621 and 632 below their current upper limits (obtained from atmospheric and solar 
experiment data). 

From the observation of (p^, if CPTV is dominant, we found that there could be large 
deviations with respect to the pure standard-oscillation case, depending on the values of 
the CPTV parameters. These deviations start at A = 0.1 (b 32 ~ 10~ 28 GeV, b 2 i ~ 10~ 26 
GeV) and reach a plateau at A = 1. There are overlaps between the different neutrino 
production models, in such a way that a measurement of certain values of (p^ could be 
satisfied by two different production models, either including CPTV or not. 

When we consider the possibility of detecting <j)^ and <p e , from which R = </>^/</> e can 
be built, we find that the regions corresponding to 1 : 2 : and 0:1:0 exhibit a similar 
behaviour than when <fi^ alone is measured, though the former is now nearly contained by 
the latter. This is not the case for 1 : : 0, where the value of R could blow up owing to a 
potentially very low value of <p e , in comparison to <f>^. 

If tau-neutrinos can also be detected, then we can use the ratio S = ^ T /0 M . When we 
combine R and S, i.e., when we assume the ability to measure <f) T with enough statistics, 
we improve the chances of discovering CPTV effects. In fact, large CPTV regions of (R, S) 
values that are well distinguished from the standard-oscillation case are obtained for the 
three neutrino production models that we have considered here. As a consequence of the 
wideness of these regions, there are many overlapping areas where a given pair (R, S) can 
be generated by any of the production models that we have explored, assuming a dominant 
CPTV. 

On the other hand, in the case of the non-observation of deviations in the flavour ratios, 
it will be possible to impose very stringent limits on the parameters related to CPTV in 
the neutrino sector, such as 621 ^ 10~ 29 GeV and 632 < 10~ 27 GeV, to be compared with 
the current limits of 10~ 23 GeV and 10~ 21 GeV, respectively. 

In order to compare CPTV with other competitive new physics scenarios, we have 
included the possibility of neutrino decay to invisible products, both in the normal and 
inverted mass hierarchies. As a result, there are additional overlaps between the CPTV 
regions and the regions accessible by neutrino decays, for certain values of the CPTV 
parameters. 

With the purpose of presenting a more realistic perspective, we have performed an 
analysis of the potential CPTV signals at a large ice Cerenkov detector such as IceCube. 
On top of the three neutrino production models, 1:2:0,0:1:0 and 1 : : 0, we have used 
two different models for the neutrino astrophysical flux, one by Waxman and Bahcall and 
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the other by Koers and Tinyakov, and compared their respective signals at the detector. 
For this analysis, we have found that there are still overlaps among the three production 
models, even more pronounced than the ones that were obtained in the theoretical plots. 
In 15 years of exposure time, a 1 km 3 detector (or, equivalently, tV e g = 15 yr km 3 ) would 
be able to distinguish between the two fluxes, in the 1:2:0 and 0:1:0 scenarios, while 
a 5 km 3 detector with the same exposure (or tV e g = 75 yr km 3 ), would provide clearer 
separation between the flux models at the la level. 

When A > 1, a separation of a few standard deviations between the CPT-conserving 
and the CPTV scenarios is possible, depending on the values of the CPTV mixing param- 
eters, both for the Waxman-Bahcall and the Koers-Tinyakov fluxes. This separation is 
more clearly visible in the 1:2:0 and 0:1:0 cases. Our main result is that, while it is in 
principle possible to detect the presence of CPTV with IceCube, it will not be possible, for 
many values of (R, S), to find which one of the production models is the actual one. The 
detection of CPTV can be aided by a more precise knowledge of the standard-oscillation 
mixing parameters and by a considerably higher event yield, brought by a larger effective 
volume or probably from a non-Cerenkov detector, but these two improvements do not 
eliminate the overlaps that exist between (R, S) regions associated to different production 
models. In the event of measuring values of R, S, or both, that fall inside an overlap region, 
then, in the absence of knowledge of the values of the CPTV mixing parameters, one is able 
to trade one production model, with a certain set of values of the CPTV parameters, by 
one of the other overlapping models, with another set of values of the parameters. These 
degeneracies could be lifted by an independent measurement of the CPTV parameters in 
the neutrino sector, in a different kind of experiment. 

Nevertheless, as a tool for detecting the presence of CPT breaking in the neutrino 
sector, if not for measuring it in detail, IceCube might be a useful one. If CPT is broken, 
and if the parameters introduced by the breaking have certain values, then IceCube could 
be able to detect the deviation from the standard-oscillation scenario after 15 years of data 
taking. 



Acknowledgments 

This work was supported by grants from the Direccion Academica de Investigacion of the 
Pontificia Universidad Catolica del Peru (projects DAI-4075 and DAI-L009 [LUCET]) and 
by a High Energy Latinamerican-European Network (HELEN) STT grant. MB acknowl- 
edges the hospitality of IFIC during the development of this work. 



References 

[1] R. J. Davis, D. S. Harmer and K. C. Hoffman, Search for neutrinos from the sun, Phys. Rev 



Lett. 20 (1968) 1205 



[2] B. Aharmim et al. [SNO Collaboration], Electron energy spectra, fluxes, and day-night 



asymmetries of B-8 solar neutrinos from the 391-day salt phase SNO data set, Phys. Rev. C 



72 (2005) 055502 [rmcl-ex/0502021 



- 24 - 



C. Arpesella et al. [Borexino Collaboration], Direct Measurement of the Be-7 Solar Neutrino 



Flux with 192 Days of Borexino Data, phys. Rev. Lett. 101 (2008) 091302 



| astro-ph/0805 . 3843 



Y. Ashie et al. [Super-Kamiokande Collaboration], Evidence for an oscillatory signature in 
atmospheric neutrino oscillation, Phys. Rev. Lett. 93 (2004) 101801 hep-ex/0404034]. 

M. Ambrosio et al. [MACRO Collaboration], Measurements of atmospheric muon neutrino 



oscillations, global analysis of the data collected with MACRO detector, Eur. Phys. J. C 36 



(2004) 323. 



K. Eguchi et al. [KamLAND Collaboration], First results from KamLAND: Evidence for 
reactor anti-neutrino disappearance, Phys. Rev. Lett. 90 (2003) 021802 |hep-ex/0212021]. 

M. H. Ahn et al. [K2K Collaboration], Measurement of Neutrino Oscillation by the K2K 
Experiment, phys. Rev. D 74 (2006) 072003| ]hep-ex/0606032 ]. 



P. Adamson et al. [MINOS Collaboration], Measurement of Neutrino Oscillations with the 
MINOS Detectors in the NuMI Beam, phys. Rev. Lett. 101 (2008) 131802| 



[hep-ex/0806.2237 



C. W. Walter, Experimental Neutrino Physics, |hep-ex/0810.3937[. 

G. L. Fogli, E. Lisi, A. Marrone and G. Scioscia, Testing violations of special and general 
relativity through the energy dependence of v T oscillations in the Super-Kamiokande 

atmospheric neutrino experiment, Phys. Rev. D 60 (1999) 053006 [hep-ph/9904248]. 



N. E. Mavromatos, Neutrinos and the phenomenology of CPT violation, [iep-ph/040200E]. 



N. E. Mavromatos, CPT Violation and Decoherence in Quantum Gravity, J. Phys. Conf. Ser. 



171 (2009) 012007 [tiep-ph/0904. 0606[. 



D. Colladay and V. A. Kostclecky, Lorentz-violating extension of the standard model, Phys 



Rev. D 58 (1998) 116002j Hhep-pti/9809521 ]. 



V. A. Kostclecky and M. Mcwes, Lorentz and CPT violation in the neutrino sector, Phys. 



Rev. D 70 (2004) 031902; Hhep-pti/0308300 ]. 



V. Dc Sabbata and M. Gasperini, Neutrino Oscillations In The Presence Of Torsion, Nuovo 



Cim. A 65 (1981) 479 



M. Bustamante, A. M. Gago and C. Pena-Garay, Extreme scenarios of new physics in the 



UHE astrophysical neutrino flavour ratios, J. Phys. Conf. Ser. 171 (2009) 012048 



|hep-ph/0906.5329 



D. Decamp et al. [ALEPH Collaboration], Determination of the Number of Light Neutrino 
Species, phys. Lett. B 231 (1989) 5jf . 

C. Amsler et al. [Particle Data Group], Review of particle physics, Phys. Lett. B 667 (2008)| 



B. Kayser, Neutrino Mass, Mixing, and Flavor Change, tiep-ph/0804. 1497. 



T. Schwetz, M. Tortola and J. W. F. Valle, Three-flavour neutrino oscillation update, New J 



Phys. 10 (2008) 113011 hep-ph/0808 . 2016|. 



[21] A. Dighe and S. Ray, CPT violation in long baseline neutrino experiments: a three flavor 
analysis, phys. Rev. D 78 (2008) 036002j ]hep-pti/0802 . 0121]]]. 



- 25 - 



J. L. Bazo, M. Bustamante, A. M. Gago and O. G. Miranda, High energy astrophysical 



neutrino flux and modified dispersion relations, Int. J. Mod. Phys. A 24 (2009) 5819 



| hep-ph/0907. 1979 



G. Barenboim and C. Quigg, Neutrino observatories can characterize cosmic sources and 
neutrino properties, \Phys. Rev. D 67 (2003) 073024| |hep-ph/0301220|| . 



J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa and T. J. Weiler, Measuring flavor ratios of 



high-energy astrophysical neutrinos, Phys. Rev. D 68 (2003) 093005 [Erratum-z&id. 72 (2005) 
019901] flhep-ph/0307025 l. 

A. Esmaili and Y. Farzan, An Analysis of Cosmic Neutrinos: Flavor Composition at Source 



and Neutrino Mixing Parameters, \Nucl. Phys. B 821 (2009) 197| Jhep-pti/0905 . 0259 . 
H. Athar, C. S. Kim and J. Lee, The intrinsic and oscillated astrophysical neutrino flavor 



ratios, Mod. Phys. Lett. A 21 (2006) 104S [hep-pti/0505017| 



P. Lipari, M. Lusignoli and D. Meloni, Flavor Composition and Energy Spectrum of 
Astrophysical Neutrinos, \Phys. Rev. D 75 (2007) 123005| ]astro-ph/0704 . 071S| ]. 



K. C. Lai, G. L. Lin and T. C. Liu, Determination of the Neutrino Flavor Ratio at the 



Astrophysical Source, Phys. Rev. D 80 (2009) 103005 [hep-ph/0905 .4003 



Z. Z. Xing and S. Zhou, Towards determination of the initial flavor composition of 



ultrahigh- energy neutrino fluxes with neutrino telescopes, Phys. Rev. D 74 (2006) 013010 



astro-ph/0603781 



V. D. Barger, S. Pakvasa, T. J. Weiler and K. Whisnant, CPT odd resonances in neutrino 



oscillations, \Phys. Rev. Lett. 85 (2000) 5055| |hep-ph/0005197 ]. 

J. P. Rachen and P. Meszaros, Photohadronic neutrinos from transients in astrophysical 



Phys. Rev. D 58 (1998) 123005| [ |astro-pri/9802280 |. 



T. Kashti and E. Waxman, Flavoring astrophysical neutrinos: Flavor ratios depend on 



energy, \Phys. Rev. Lett. 95 (2005) 18110l| [ |astro-ph/0507599 | 



Z. z. Xing, Neutrino telescopes as a probe of broken mu tau symmetry, Phys. Rev. D 74 
(2006) 0130"09| jriep-ph/0605219 1. 



M. Maltoni and W. Winter, Testing neutrino oscillations plus decay with neutrino telescopes, 
J. High Energy Phys. 0807 (2008) 064| [ |hep-ph/0803 . 205(fl . 



J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa and T. J. Weiler, Decay of high-energy 



astrophysical neutrinos, Phys. Rev. Lett. 90 (2003) 181301 [hep-ph/0211305|. 



J. F. Beacom and N. F. Bell, Do solar neutrinos decay? Phys. Rev. D 65 (2002) 113009 



[hep-ph/0204111]. 

A. Mirizzi, D. Montanino and P. D. Serpico, Revisiting cosmological bounds on radiative 



neutrino lifetime, Phys. Rev. D 76 (2007) 053007 hep-ph/0705 . 4667 



M. Blennow and D. Meloni, Non-standard interaction effects on astrophysical neutrino fluxes, 
Phys. Rev. D 80 (2009) 065009| ftiep-ph/0901 . 2110| . 



S. Ando, M. Kamionkowski and I. Mocioiu, Neutrino Oscillations, Lorentz/CPT Violation, 



and Dark Energy, \Phys. Rev. D 80 (2009) 123522| flriep-ph/0910 . 4391 



Z. z. Xing and S. Zhou, Cosmic Neutrino Flavor Democracy and Unitarity Violation at 
Neutrino Telescopes, \Phys. Lett. B 666 (2008) 166| friep-ph/0804 . 3512 |. 



- 26 - 



T. DeYoung, Neutrino Astronomy with IceCube, Mod. Phys. Lett. A 24 (2009) 1543 



| astro-ph/0906 . 4530 | . 

F. Halzen, IceCube: The Rationale for Kilometer- Scale Neutrino Detectors, 



[astro-ph/0910.0436]. 

L. A. Anchordoqui, H. Goldberg, M. C. Gonzalez-Garcia, F. Halzen, D. Hooper, S. Sarkar 



and T. J. Weiler, Probing Planck scale physics with IceCube, Phys. Rev. D 72 (2005) 065019 



[ hep-ph/0506168 [. 

R. Gandhi, C. Quigg, M. H. Reno and I. Sarcevic, Neutrino interactions at ultrahigh- energies, 



Phys. Rev. D 58 (1998) 093009| ^iep-ph/9807264]l 



R. Gandhi, C. Quigg, M. H. Reno and I. Sarcevic, Ultrahigh- energy neutrino interactions, 



Astropart. Phys. 5 (1996) 81 [hep-ph/9512364| 



E. Waxman, Astrophysical sources of high energy neutrinos, Nucl. Phys. 118 (Proc. Suppl 



(2003) 353 astro-ph/0211358| 



J. K. Becker and P. L. Bicrmann, Neutrinos from active black holes, sources of ultra high 
energy cosmic rays, \Astropart. Phys. 31 (2009) 138| Jastro-pti/0805 . 1498[ . 



H. B. J. Koers and P. Tinyakov, Relation between the neutrino flux from Centaurus A and 



the associated diffuse neutrino flux, \Phys. Rev. D 78 (2008) 083009| flastro-ph/0802 . 2403 



A. Achterberg et al. [IceCube Collaboration], Multi-year search for a diffuse flux of muon 



neutrinos with AMANDA-II, Phys. Rev. D 76 (2007) 042008| [Erratum- i&id 77 (2008) 
089904] ]astro-ph/0705. 13151 ~ 



-27- 



